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Abstract 

The MiniBooNE and LSND experiments are compatible with each other when two 
sterile neutrinos are added to the three active ones. In this case there are eight 



possible mass orderings. In two of them both sterile neutrinos are heavier than the 



<^> three active ones. In the next two scenarios both sterile neutrinos are lighter than 

the three active ones. The remaining four scenarios have one sterile neutrino heavier 
and another lighter than the three active ones. We analyze all scenarios with respect 
to their predictions for mass-related observables. These are the sum of neutrino 
^ ■ masses as constrained by cosmological observations, the kinematic mass parameter as 

measurable in the KATRIN experiment, and the effective mass governing neutrinoless 
double beta decay. It is investigated how these non-oscillation probes can distinguish 
between the eight scenarios. Six of the eight possible mass orderings predict positive 
signals in the KATRIN and future neutrinoless double beta decay experiments. We 
also remark on scenarios with three sterile neutrinos. In addition we make some 
comments on the possibility of using decays of high energy astrophysical neutrinos 
to discriminate between the mass orderings in presence of two sterile neutrinos. 
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1 Introduction 



The long awaited results of the MiniBooNE experiment pQ showed that scenarios in which 
one sterile neutrino is added to the three active ones are incompatible with data. Such 
schemes were motivated by the results from the LSND experiment [2] , which observed flavor 
transitions interpreted as <-> v e neutrino oscillations. A number of authors [3j HJ El [7] 
investigated the implications of such schemes. In particular, it was realized [8] that so- 
called 2+2 scenarios (two pairs of neutrinos close in mass separated by a large gap) are 
ruled out and that only 3+1 scenarios (three mostly active neutrinos separated by a large 
gap from the mostly sterile one) are allowed, though only small part of the parameter 
space survived. The MiniBooNE results ruled out even this part [9] by excluding the 
LSND parameter space at the 98 % C.L. []]. 

Allowing one more sterile neutrino to enter the stage improves the compatibility of LSND 
with other experiments [TQj, |TT] and in particular renders the MiniBooNE and LSND expe- 
riments compatible [9]. Only comparably few models for neutrino schemes with two extra 
sterile neutrinos have been constructed [12], and the potentially rich phenomenology of 
such scenarios is hardly investigated [13], [H] . 

With two sterile neutrinos added to the usual three, one has eight possible mass orderings, 
which should be compared with the two schemes (normal and inverted ordering) in case 
of "only" three active neutrinos. We study in this paper the predictions of the eight cases 
for mass-related observables. We investigate the sum of neutrino masses as constrained by 
cosmological observations, the kinematic mass parameter as measurable in the KATRIN 
experiment (partly analyzed also in Ref. [ID]), and the effective mass controlling neutri- 
noless double beta decajlE We also investigate how and if mass-related observables can 
contribute to distinguish the possibilities. The mass patterns in the 3+2 scheme can be 
classified in three main classes: 

• two 2+3 scenarios: the two sterile neutrinos are heavier than the three active ones; 

• two 3+2 scenarios: the two sterile neutrinos are lighter than the three active ones; 

• four 1+3+1 scenarios: one sterile neutrino is heavier than the three active ones which 
in turn are heavier than the second sterile neutrino. 

The paper is build up as follows: first we summarize the required formalism in Section [2] 
before outlining in Section [3] the eight possible mass orderings for scenarios with two sterile 
neutrinos. In Section H] we study the mass-related observables for the two scenarios which 
have the sterile neutrinos heavier than the active ones, while in Section [5] the two scenarios 
in which the sterile neutrinos are lighter than the active neutrinos are analyzed. Section 
|6] contains the four 1+3+1 scenarios and a short discussion on scenarios with three sterile 
neutrinos is delegated to Appendix |A] In Appendix [B] we discuss another interesting possi- 
bility to distinguish between the different scenarios outlined above at neutrino telescopes, 

1 For related analyzes in different sterile neutrino scenarios, see, e.g., Ref. [IB] . 
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allowing high energy astrophysical neutrinos to decay. Finally, in Section [7] we discuss and 
summarize our findings. 
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2 Formalism 

2.1 Neutrino Mixing 

Neutrino mixing is described by the leptonic mixing, or Pontecorvo-Maki-Nakagawa-Sakata 
(PMNS) matrix U: 



(1) 



It links the mass eigenstates 1^1,2,3,4,5 having masses 7711,2,3,4,5 with the three active flavor 
states v e ,fj,,T and the two new sterile states v Sl and v S2 . In what follows we denote the 
mostly active neutrinos with ^1,2,3- Regardless of their ordering (normal or inverted) we 
have 

\U e i\ 2 ~ cos 2 6» Q , \U e2 \ 2 ~ sin 2 6 Q and |/7 e3 | 2 ~ sin 2 6» C hooz , (2) 

where 9 & is the mixing angle for solar and KamLAND neutrinos and #chooz the mixing 
angle for short baseline reactor neutrinos. We further have the mass-squared differences 
governing solar and atmospheric neutrino oscillations. In the following, we will use the 
following best-fit and 3a ranges [T6] : 

sin 2 # = 0.30i{j:{j| with Am l = (8.oll:8) • eV 2 , 



(3) 

sin 2 flcHooz = 0.00±8:8 4 with Ami = (2.6±g;|) • 10" 3 eV 2 



The related typical mass scales are therefore yAm Q ~ 0.009 eV and yAmj^ ~ 0.05 eV, 
respectively. In what regards the two additional sterile neutrinos, the analysis in Ref. [9] 
resulted in the following best-fit values^ 

Am 2 ! = 6.49l^g eV 2 with \U e5 \ = 0.12 , 

(4) 

Am s 2 2 = 0.89±gi eV 2 with \U eA \ = 0.11 . 

In what follows we will give for our observables explicit numerical values obtained with not 
only these best-fit values, but also for another typical illustrative point in the parameter 
space: 

Am 2 ! = 1.9018:98 e V 2 with \U e5 \ = 0.12 , 
Am 2 s2 = 0.9018:20 eV 2 with \U e4 \ = 0.11 . 



2 We consider only the analysis which uses the MiniBooNE results above reconstructed neutrino energies 
of 475 MeV, because results from the lowest energy bin are not well understood pQ. 
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The main feature of this point is of course the smaller overall neutrino mass it implies. It 
corresponds approximately to the center of another, isolated region allowed at 90 % C.L. of 
Figure 6 from Ref. [9]. The two central points from Eqs. (j4j) and (111) are quite typical for 
the situation in the presence of two sterile neutrinos and we will make frequent use of them. 
With the values of the two new mass-squared differences we can estimate typical neutrino 
mass scales, which we will encounter frequently in the following: JKm^ ~ 2.55 (1.38) 

eV, yjAm% 2 ~ 0.94 (0.95) eV, ^Am 2 sl - Am 2 s2 ~ 2.37 (1.00) eV, and ^Am 2 sl + Am 2 2 ~ 
2.72 (1.67) eV. 

In Ref. [9] the allowed ranges of the mixing matrix elements \U e i\ and \U e 5\ are not given. 
However, we will see in our discussions that especially for neutrinoless double beta decay 
it is sometimes important to analyze the impact of varying these parameters. In absence 
of any information we have varied \U e ^\ and \Ues\ by 50% around the best-fit points in 
Eqs. (jl]) and ([5]). Thus we consider the following ranges for the parameters \U e ^ and \U e <~\: 

|C/ e5 | = 0.12+° : °6 and |[/ e4 | = 0.11±°;°| . (6) 
for both the best-fit point and the second illustrative point considered in Eq. (jSJ). 



2.2 Neutrino Masses 

As neutrino oscillations are sensitive only to mass-squared differences, the neutrino mass 
scale is not known, but only limited from above by different experiments and observations. 
Typically, the mass scale is inversely proportional to the scale of the mechanism which 
is responsible for neutrino mass. Therefore, knowing the mass is a very important step 
towards the understanding of neutrino physics. 
Mass-related observables are the sum of neutrino masses 

E = E™- (7) 

i 

which can be inferred from cosmological observations. Typical limits are smaller than 
about 1 eV [T71 (TSJ dHl ED], but they depend on the used data sets, the number of neutrino 
species and how the mass is distributed among the different neutrinos. We will discuss 
this in more detail in Section [71 One also has the kinematic neutrino mass parameter 
measurable in /3-decay experiments 

m p = p^Wrn* . (8) 

This quantity is measured when the electron energy interval around the endpoint of the 
investigated beta decay is much larger than the rrij, otherwise corrections to this formula are 
required [21] . This condition is fulfilled for the values of the neutrino masses we choose^. 

3 In principle, the analysis of Ref. [9] allows sterile neutrino mass values of Am 2 > 10 eV 2 , which are 
indeed close to the energy interval used by the upcoming KATRIN experiment f22j. These values are 
however in very strong conflict with all mass-related observables and we therefore omit them. For the 
values used we estimate corrections to Eq. ([8} to be at most of order 10%. 
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NH 




yj sin 2 6 Q Ami + sin2 ^chooz Am\ 


sin 2 6 Q ^Am 2 Q + sin 2 C hooz ^Aml e^ 2 


IH 


V Am A 




^Am^ yjl- sin 2 26> sin 2 a 2 /2 


QD 


3 rriQ 


m 


m ijl — sin 2 26> sin 2 02/2 



Table 1: Extreme limits of 3-flavor scenarios and the resulting mass- related observables. 
We have defined 032 = 0:3 — 02. 



The current limit on 7713 is 2.2 eV at 95% CL. [23 
magnitude is expected by the KATRIN experiment 
mass in neutrinoless double beta decay: 



and improvement by one order of 
[22] . Finally, we have the effective 



(to) 



in, 



E 



1 2 e iai 



m, 



(9) 



Here 02,3,4,5 are the four possible and unknown^ Majorana phases (we can choose 01 =0). 
Applying nuclear matrix element uncertainties on the current 90% CL. limits on lifetimes 
[2H [25] gives limits on the effective mass in the range of 1 eV, and sizable improvement is 
expected also in this field [26J. 

Let us summarize for the sake of comparison the situation in 3-flavor scenarios (see for 
instance [27, 28J). We have three extreme cases of the mass ordering, the normal hierarchy 
(NH, toJ ~ Am\ 3> to 2 ; — Ato ^> to 2 ), the inverted hierarchy (IH, to 2 ~ to 2 ~ Am\ ^> 
to|) and quasi-degenerate neutrinos (QD, to| ~ to 2 — to 2 = TOq ^> Ato^Atoq). Table 
[T] shows the results for the mass-related observables. Obviously, to^ is unobservably low 
for NH and IH, while (to) and E are for NH. It may be possible to probe the inverted 
hierarchy regime through future cosmo logical observations |17j . 



Certain mass orderings to be discussed in the following will have problems with some of 
the three observables E, mp or (to). Does this mean that they are ruled out? Not neces- 
sarily, because cosmological neutrino mass (and number [3D]) limits can in principle 
be evaded, or relaxed by a factor of a few, by means of unknown neutrino interactions or 
other cosmological features. Furthermore, confronting a mass ordering with the limit on 
the effective mass makes only sense when neutrinos are Majorana particles, which however 
is a very well justified assumption. Only the kinematic parameter mp does not suffer from 
any underlying model assumption and provides an unambiguous test. We will leave aside 
discussions of the validity of the different limits in particular on E. Our aim is simply to 
study the predictions for the observables for all eight possible mass orderings, which we 
will outline in the next Section. 

4 In fact, the analysis of neutrino scenarios with two sterile neutrinos gives some constraint on a CP 
phase [9] (see also [H]). Being a "Dirac-phase" , however, it does not appear in (m). 
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3 Eight Possible Mass Orderings in Neutrino Scenar- 
ios with two sterile Neutrinos 



As mentioned in the Introduction, analyzes of the LSND and MiniBooNE, as well as 
various other experiments, give a consistent picture only if two additional independent 
mass-squared differences Am^ > Am s 2 2 are present. We assume here that the difference 
between the two, Am 2 ss = Am^ — Am^, is much larger than Am\, an assumption used 
in the 3+2 analyzes in [9] from where we take the values of the additional parameters. 
Our convention is the following: we call the masses of the three predominantly active 
neutrinos 777.3, m 2 and mi. The mixing among them is responsible for the solar, atmospheric 
and short baseline reactor neutrino oscillation results. If they are normally ordered, then 
m 3 > m 2 > rrii with 

m\ — m\ = Am\ and m\ — m\ = Am^ . (10) 

They can also be inversely ordered, in which case 777,2 > rri\ > m$ and 

m\ — m\ = Am\ and m\ — ml = Am 2 & . (11) 

We have to add two predominately sterile neutrino states, whose masses we denote by 
7774 and 777.5. They can either be heavier or lighter than the three active neutrinos (in 
what follows, we will omit for simplicity the word "predominately" or "mostly" in front 
of "sterile neutrinos" and "active neutrinos"). Without loss of generality, we can choose 
that 7775 is either the largest or the smallest mass and associate it always with Am^, while 

7774 is always associated with Am s 2 2 . In this case, we do not have to rename the matrix 
elements C/ e 4 and [/ e 5, which quantify the mixing of the two additional neutrinos with the 
electron neutrino. In general, one could choose the labels of the masses such that always 

7775 > 7774 > 7773 > r?7 2 > mi holds. In this case, however, the mixing matrix elements 
would be different for each of the possible mass orderings. With our convention, the values 
of \U e i\ are fixed by Eqs. (T2]) and (jlj [5]) and do not have to be relabeled. 

Let us first discuss the case of both sterile neutrinos being either heavier or lighter than the 
three active ones: the largest of the two independent new mass-squared differences, Am s 2 j, is 
then always the largest possible mass-squared difference. If the two sterile states are above 
the three active ones ("2+3 scenarios"), then Am 2 2 is the mass-squared difference between 
the lightest sterile state (which is also the second heaviest state) and the lightest available 
state (which is active). If the two sterile states are below the three active ones ("3+2 
scenarios"), then Am 2 2 is the mass-squared difference between the heaviest state (which is 
active) and the heaviest sterile state (which is the second lightest state). The four possible 
schemes are shown in Figs. [1] and [2j The names of the schemes are defined as follows: 
depending on whether the active neutrinos are normally or inversely ordered, the scheme 
has the capital letter "N" or "I" in its name. Depending on whether the two sterile neutrino 
masses 7774 and 7775 are lighter or heavier than the active ones, the capital letters "SS" appear 
after or in front of this capital letter. For instance, if 7775 > 7774 > 7773 > 7772 > m±, then we 
call the scenario SSN, while for 7772 > mi > m 3 > m 4 > m 5 we call it ISS. 
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The last possible class of mass orderings is when one sterile neutrino is heavier then the 
three active ones, which in turn are heavier than the second sterile neutrino ("1+3+1 sce- 
narios"). The heaviest neutrino can either be separated by Am^ or Am s 2 2 from the active 
neutrinos. If it is separated by Am s 2 j (Am s 2 2 ) then we call the scenario SNSa or SISa (SNSb 
or SISb). The possibilities are shown in Fig. [3J We note here that the fit of Ref. [5J, and also 
the analyzes of Refs. p^Ol [TTj . do strictly speaking not apply to these schemes. The reason 
is that in the oscillation probabilities for — > v e transitions there are not only terms pro- 
portional to sin 2 Aro 2 j and to sin 2 Am 2 2 but also an interference term proportional 
to cos(Am 2 + 5), where 6 is a CP phase (which does not appear in survival probabil- 
ities). Refs. [HI QUI ELI] make the implicit assumption |Am 2 | = lAm 2 ^ — |Am 2 2 |, which is 
not fulfilled for the 1+3+1 scenarios. In lack of any detailed fit of the data within these 
schemes we will assume for simplicity that the mass-squared differences are the same. The 
values of the mass-squared differences resulting from fits taking into account the 1+3+1 
case will remain of course in the eV range. 

Summarizing we end up with eight different schemes^]. We can already at the present stage 
make some general statements. First of all, the sum of neutrino masses depends basically 
only on the new mass-squared differences and typical values will be 



E > sjAm 2 sl + VAm s 2 2 or E > 3 ^Am 2 8l or E > 3 ^Am 2 s2 , (12) 

depending on the details of the mass ordering. In general, all of them are expected to 
face serious problems with cosmology, and the smaller the mass-squared differences Am s 2 j 
and Am 2 2 are, the smaller E. Another point worth mentioning is that the effective mass 
governing neutrinoless double beta decay can be written as 

(m) = (m) 3 + (m) £ 
where (m) 3 = cos 2 9 & ni\ + sin 2 9 Q m 2 e ta2 + sin 2 #chooz ^3 e* a3 (13) 
and (m) st = | U e4 | 2 m 4 e iai + \U e5 \ 2 m 5 e iOLh . 

Obviously, |(^) 3 | is an effective mass similar to the one analyzed in the usual three-flavor 
situation [27], cf. Table [U The quantity \(m) s \ is the contribution from the two sterile 
states. We will encounter in the following all cases: dominance of the sterile contribution, 
dominance of the active contribution, and equal-sized contributions, leading potentially to 
complete cancellation. The same cases are also present for the kinematic neutrino mass 
171/3, where however no cancellation is possible as it is given by an incoherent sum. 

We will discuss now the three mass-related observables for the eight possible mass orderings. 
We give approximate analytic expressions for these observables in the limit of vanishing 
smallest mass and sin 2 6*cHooz- We use these expressions to give illustrative numerical 
values in each case for the best-fit values of the oscillation parameters (or for the second 



5 In scenarios with three sterile neutrinos one would have 16 possible mass orderings, see Appendix 1X1 
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illustrative central point from Eq. ((Sj)). We will also plot the observables as a function of 
the smallest neutrino mass for the central points as well as by varying the parameters in 
their corresponding allowed ranges from Eqs. (El HI [5]). 



4 Sterile Neutrinos heavier than active Neutrinos: 2+3 
Scenarios 

4.1 Scheme SSN 

In this scheme, 777,5 > m 4 > m 3 > m 2 > m \- The three lowest states account for the solar 
and atmospheric neutrino mass-squared differences according to Eq. ffTUl) . We have 

Schematically, this scheme is shown in Fig. [TJ We can express the individual masses in 
terms of the smallest mass mi and the independent mass-squared differences: 



m2 = J Amg + m 2 , m.3 = J Am\ + Am 2 & + mf , 



m 4 = J Am 2 2 + m{ , 777,5 = 1/ Am s i + m{ 



(15) 



The typical masses are therefore 7772 — \j Am & ~ 0.01 eV, 7773 ~ y Am A ~ 0.05 eV, 

m 4 ~ y/Am% 2 ~ 0.94 (0.95) eV and m 5 ~ ^Am^ ~ 2.55 (1.38) eV. The three lightest 
neutrinos have the same values as in a normal hierarchical 3-flavor framework. The upper 
left plot in Fig. H] shows the individual masses as a function of the smallest mass in this 
picture. The limit where all five neutrinos are quasi-degenerate comes when the smallest 
mass is beyond 1 eV. 

If this scheme is realized then the sum of neutrino masses which is constrained from cos- 
mology is given as 



-ssn ^ lAm% + yAmi + \J Am 2 sl + \j Am 2 s2 

(16) 



Am^ + ^jAm 2 s2 ~ 3.49 (2.32) eV . 

Neglecting mi is a good approximation as long as mi ^0.1 eV. For such values, E lies 
roughly between 3.3 and 3.8 (or 1.9 and 2.6) eV, but it can reach unrealistically large 
values of 10 eV for a smallest mass in the eV range. 

The kinematic mass is also mainly given by the sterile neutrino contribution: 



m 



SSN 

P 



sin 2 6 & Am| + sin 2 # C hooz Am| + \U ei \ 2 Am 2 2 + \U e5 \ 2 Am 2 sl 



U e 4 2 Am 2 2 + \U e5 \ 2 Am 2 sl ~ 0.32 (0.20) eV . 



(17) 



Both S and mp are shown in Fig. [5], where the solid lines are these quantities at the best-fit 
values of Eq. (jlj) whereas the bands are obtained by varying the masses and the mixing 
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angles within their allowed ranges. The KATRIN experiment, having a sensitivity of 0.3 
eV, will find a positive signal if m\ ~?t 0.3 eV. For smaller values, however, mp can lie below 
0.3 eV. If mp larger than 0.5 eV is found then this scenario is ruled out, unless mi ^ 0.2 
eV. No qualitatively different features are found for the second point from Eq. (jHJ). 
Finally, neutrinoless double beta decay should be triggered by an effective mass given by 



/ \SSN 

(m) ~ 



sin 2 6 & J Am 2 + sin 2 # C hooz J Am\ e^ 3 "^ 



+ \U e4 \ 2 JAm 2 2 e i( - ai ~ a ^ + \U e5 \ 2 J Am 2 sl e *(«*-«a 



U ei \ 2 jAm 2 2 + \U e5 \ 2 jAm 2 sl e^"^) 



(0.025 0.048) eV 



where the two sterile neutrinos provide the leading contribution. In case of the second 
typical point from Eq. (J^j) the two additional neutrinos give a leading contribution between 
0.008 and 0.031 eV. The upper left panels of Figs. M and show the effective mass as a 
function of the smallest mass in this scenario. The shaded region inside is drawn for the 
best-fit values of mass and mixing parameters and varying the Majorana phases between 
and 27r, while for the outer shaded regions we vary these parameters also in their permissible 
ranges. The figures show that neglecting the smallest mass is a good approximation as long 
as mi ~ 0.01 eV. There is for both central points a cancellation regime (for mi between 
0.02 and 0.1 eV) in which the effective mass vanishes or becomes unobservably small. If 
we use the ranges around the two central points of the sterile neutrino parameters, then 
the two terms can cancel when the conditions \U e ±/U t 



(:5 



Am 2 i/Am 2 2 and a 5 — a 4 



7T 



are fulfilled. For the smallest neutrino mass above 0.1 eV the effective mass cannot vanish 
due to the non-maximal solar neutrino mixing angle. For smaller values of mi the scenario 
is ruled out if (m) is found to be larger than 0.1 eV. 



4.2 Scheme SSI 

In this scenario, schematically shown in Fig. [TJ it holds m 5 > m 4 > m 2 > mi > m 3 , i.e., 
the two heavy sterile neutrinos are heavier than the three light neutrinos which enjoy an 
inverted hierarchy. Consequently, Eq. (Till holds. In addition, we have 



Am 2 ! 



m 5 



m\ and Am 2 2 



m. 



m 



3 i 



and the masses in terms of the smallest mass are 



mi = y Am\ + m 2 , m 2 = J Am\ + Am 2 Q + m 
m 4 = \J Am 2 2 + m 2 , m 5 = \J Am s \ + m 2 . 

mi 



3 • 



Am 2 > 



(19) 



(20) 



0.05 eV, 



Neglecting the smallest mass we find typical values of m 2 

0.94 (0.95) eV and m 5 ~ jAm 2 sl ~ 2.55 (1.38) eV. The active neutrinos 



m 4 



'Am 



s2 



behave according to an inverted hierarchy in a 3-generation framework. As a function of 
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the smallest neutrino mass 7773 the other masses are shown in Fig. HI One finds that S SSI 
is identical to S in Eq. (fT6|) (Am has to be replaced with Am\, which does not make 
a notable difference). The kinematic mass is also basically identical to the one in scenario 
SSN, which is given in Eq. (1171) . In the effective mass the situation is different^, because 
the light neutrinos obey an inverted hierarchy and therefore a contribution of the same 
order of magnitude as the sterile ones: 



/ \SSI 

(m) ~ 



sin 2 # J Am 2 A + c °s 2 V Am A e^ 2 



+ \Uei\ 2 J&m 2 s2 e iQ4 + \U e5 \ 2 jAm 2 sl e ia5 



(21) 



The absolute value of the first two terms is (see Table OQ) yj Am\ sjl - sin 2 28 Q sin 2 a 2 /2 
(between 0.020 and 0.051 eV) while for the best-fit values from Eq. PJ the absolute value 
of the last two terms is between 0.025 and 0.048 eV, see Eq. (|18|) . Hence, the effective mass 
can vanish completely in this scheme even for the best-fit values and a vanishing smallest 
neutrino mass. This is borne out by the dark shaded regions in the upper right panels of 
Figs. [6] and [7] for the best-fit points of Eqs. (J4]) and ([5]). The effective mass can be as large 
as 0.1 eV (or 0.08 eV for the second typical point from Eq. ([5])) in the small mi (^ 0.1 eV) 
regime. Finding a larger (m) will rule out scenario SSI. 

We remark here that only in the schemes SSN and SSI the magnitudes of ?7 e 4 and U e s 
are important for the predictions of m@ and (m). However, our statements regarding the 
possible exclusion of scenarios SSN and SSI with future measurements is rather insensitive 
to the precise values of U ei and U e5 . 



5 Sterile Neutrinos lighter than active Neutrinos: 3+2 
Scenarios 

5.1 Scheme NSS 

In this scheme, m 3 > m 2 > vn,\ > 7774 > 777,5, i- e -; the three normally ordered active 
neutrinos are heavier than the two sterile neutrinos, see Fig. [2j Apart from Eq. (TTOT) it 
holds that 

Am s 2 j = m\ — m\ and Am 2 2 = m\ — m\ , (22) 
and the masses in terms of the smallest mass and the mass-squared differences are 

7774 = \j Am s 2 j — Am\ 2 + m\ , mi = sj Am^ — Am\ — Am 2 & + m\ , 

I / 

m2 = J Am s 2 j — Am\ + m 2 , 7713 = J Am^ + m\ . 

6 There is another non-oscillation probe which can distinguish SSN and SSI, namely the decay of astro- 
physical high energy neutrinos, treated in Appendix iBl 
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We therefore have for a negligible smallest mass three quasi-degenerate neutrinos and 
another quite massive state. Their values are mi ~ m 2 ~ m 3 ~ J Am s 2 j ~ 2.55 (1.38) eV 

and 1714 — J ^ m li — Am 2 2 — 2-37 (1.00) eV. As a function of the smallest mass ms, they 
are given in Fig. HI 

Let us neglect the smallest neutrino mass and insert the best-fit values. In this case, the 
sum of masses as constrained by cosmological observations is given by 

S NSS ~ ^Am 2 sl - Am s 2 2 + 3^ Am 2 sl ~ 10.0 (5.14) eV . (24) 

The kinematic mass is given by 

mjf s ~ ^ (1 - \U e5 \ 2 ) ~ 2.55 (1.38) eV . (25) 

With \U e 5\ 2 being very small there is basically no dependence on the mixing parameters 
\U e i\ 2 because of the quasi-degenerateness of the three leading neutrinos. Only for 
reaching eV values one observes deviations from the last two equations in Fig. [51 The cos- 
mological observable £ is quite large, namely between 9 and 11 (3 and 6) eV. Interestingly, 

mp ~ yAm^ lies always above 1 eV, and is even above the current bound of 2.3 eV for the 
best-fit point. Hence, this scenario constraints Ara s 2 j to lie below ~ 5.3 eV 2 . In general, 
if this scenario is realized, KATRIN will definitely observe a positive signal the absence of 
which can rule out this scenario. 

Finally, for neutrinoless double beta decay we have (in the limit \U e ^\ — > 0) 

(m) NSS ~ ^ y/l - sin 2 29 Q sin 2 a 2 /2 ~ (1.02 4- 2.55) (0.55 -=- 1.38) eV , (26) 

which is dominated by the three active quasi-degenerate neutrinos and cannot vanish due to 
the non-maximality of solar neutrino mixing as is reflected in the lower left panel of Fig. [6l 
The effective mass ranges from JAm s \ cos 29 Q to \J Am^ . In general, the effective mass 
is sizable in scenario NSS and already a part of it around the best-fit point is disfavored 
by the current upper limit of 1 eV. In fact, improving the limit on the effective mass below 
0.2 eV rules out this scheme if neutrinos are Majorana particles. 

It is possible to set limits on Am^ and in particular on 9 & and the Majorana phase a 2 
demanding (m) to lie within a specific limit. Using the 3a ranges from Eq. ([3]) and the 
ranges around the best-fit point from Eq. (jlj) one can investigate what values are allowed. 
The result for ct 2 and sin 2 9 Q can be seen in Fig. [HI We took for (m) the current limit of 
1 eV and a future limit of 0.5 eV. Taking the second central point for the sterile neutrino 
parameters gives hardly any constraint for a limit of 1 eV, but for (m) < 0.5 eV the plot 
looks similar to the 1 eV plot of the best-fit point. 

5.2 Scheme ISS 

In this scheme (see Fig. [5]) it holds m 2 > mi > m 3 > m 4 > m^, i.e., the three inversely 
ordered active neutrinos are heavier than the two sterile neutrinos. Apart from Eq. (TTTj) 
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we have 



Am^j = m\ — m\ and Am^ = m\ — m\ , (27) 
and the masses in terms of the smallest mass are 



WI4 = y Am sl — Am^ 2 + mi, , m% = y ^ m ii ~ Am A — Am Q + m 5 , 



mi = J Am^ — Arfig + m\ , m 2 = J Am^ + m| . 



(28) 



The results are basically identical to scenario NSS. We therefore do not give the expressions. 
Distinguishing between scenarios ISS and NSS could for instance be done via matter effects 
in oscillation experiments or in supernovae |13j . 

6 One heavy and one light sterile Neutrino: 1+3+1 
Scenarios 

We will discuss now the mass-related observables when the three active neutrinos are 
'"sandwiched" between the sterile ones. Recall that the fit from [9] does not apply in this 
case. We note that apart from the mass-squared differences the mixing matrix elements £7 e 4 
and £7 e 5 might also be different. However, their values do hardly influence the predictions. 
We will insert for the rest of this Section the numerical values from Eqs. (jH [5]) for the 
sterile neutrino parameters, but will indicate that there might be differences by replacing 
in the expressions Am^ — > Am s 2 j and Am^ — > Am^. 

6.1 Scheme SNSa 

In this scheme, m 5 > m 3 > m 2 > mi > m 4 . Apart from Eq. (flQl) we have 

Am s 2 j = m\ — ml and Am^ = m\ — m\ , (29) 
see Fig. [31 We can express the individual masses as 



mi = \l Am^ 2 + ml , m 2 = J Am s 2 2 + Am| + m| , 



m 3 = J Amg 2 + Amg + Am\ + ml, m 5 = J Am s \ + Am^ + ml 



(30) 



We therefore have three quasi-degenerate active neutrino masses of order y Am 2 2 ~ 0.94 (0.95) 



eV and one very heavy mass around y Am^ + Am^ 2 ~ 2.71 (1.67) eV. A plot of the m^ as 
a function of m^ is given in Fig. HI Neglecting the smallest mass is correct as long as it is 
below 0.5 eV. We can estimate that 



E SNSa _ sJAm 2 s2 + jAmli + Am s 2 2 ~ 5.55 (4.52) eV . (31) 
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The magnitude of E is quite sizable, and can be between 5 and 6 (4 and 5) eV, if one varies 
the mass-squared differences in their allowed ranges. 

The sterile masses provide the leading contribution not only in this observable, but also in 
the kinematic mass: 



m SNSa _ ^Am 2 s2 + \U e5 \ 2 (Am 2 sl + Am s 2 2 ) ~ 1.00 (0.97) eV . (32) 

As in scenarios NSS and ISS, mp ~ Am 2 2 is always above the KATRIN sensitivity of 0.3 
eV, therefore a signal in this experiment corresponding to at least 1 eV should be observed 
if this scheme is realized. 

Finally, neutrinoless double beta decay should be triggered by an effective mass given by 



(to) 



cos 2 6 Q \/Am 2 2 + sin 2 6» \J Am 2 s2 e ia2 

(33) 



+ \U e5 \ 2 jAm 2 sl + Am 2 2 e ia5 ~ jAm 2 2 Jl- sin 2 26 Q sin 2 a 2 /2 



where the third term can be neglected. Hence, (m) SNSa is dominated by the active neutrinos 



and lies between y Am 2 2 and y Am 2 2 cos 29 Q , which is roughly 0.9 and 0.38 eV, respectively. 
The effective mass for this scenario is plotted in the lower right panel of Fig. O for the best- 
fit values as well as by varying all parameters within their allowed range. The situation for 
(m) is unfortunately similar to scenarios NSS and ISS, even though here the overall mass 



scale is yAm 2 2 , while it was yAm^ in the previous cases. The problem is of course the 
allowed range of the mass-squared differences and the unknown Majorana phase. Only if 
the condition 



Am 2 lC os2fl >C \JAm 2 s2 (34) 

is fulfilled, then we can distinguish scenarios NSS/ISS and SNSa/SISa via neutrinoless 
double beta decay. In Eq. (|34|) we have included a factor £ > 1, which takes into account 
the nuclear matrix element uncertainty, a necessity when one tries to distinguish different 
mass orderings via neutrinoless double beta decay [28]. For the best-fit values and £ = 1 
indeed Eq. (I34"j) is fulfilled, but already for the second central point one cannot distinguish 
the schemes anymore as can also be seen from the lower panels of Fig. [7J 
Anyway, if neutrinos are Majorana particles, then we can rule out scenario SNSa if (m) < 
0.1 eV. One can generate plots as shown in Fig. [8] in order to obtain constraints on the 
parameters 9 Q and sin 2 a 2 from experimental information about (to). This requires limits 



which are stronger by a factor y Am^/Am^ than the limits used to generate Fig. [8J 
6.2 Scheme SISa 

In this scheme we have m 5 > m<i > mi > m 3 > 777,4, i-e., the sterile neutrinos are above 
and below three inversely ordered active neutrinos, see Fig. One finds 

Am^ = m\ — m 2 and Am 2 2 = m 2 — m\ , (35) 
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and can express the individual masses as 



777,3 = V ^ m s2 + m 4 J m l = V ^ m s2 + + 7?74 , 

(36) 



m 2 = y Ai7ig 2 + Am A + Anig + 7774 , 7775 = y Am^ + Att7 2 2 + m\ . 

We do not give the expressions for £, mp or (m), because the results are indistinguishable 
from scenario SNSa. Again, matter effects in neutrino oscillation experiments could be 
used to distinguishing the scenarios. 

6.3 Scheme SNSb 

In this scheme, 7774 > 7773 > 7772 > mi > 7775. Apart from Eq. (flOl) we have 

Ajftg! = m 2 — mg and A777g 2 = 7774 — m 2 , (37) 
see Fig. [3l The individual masses are 



777! = y Amfj + 7775 , 777,2 = V A)H s j + A777q + 777 5 , 



7773 = y Am^i + Ami + A?77 A + 7775 , 7774 = J Am^i + Amj 2 + 777 



(3* 



The mass-related observables are obtained from the formulae for scenario SNSa from Sec- 
tion [6J] with the exchange Am^ <-> Am 2 2 . Hence, 



^SNSb 



,SNSb 



~ 3JAm 2 sl + jAm 2 sl + Am 2 2 ~ 10.36 (5.81) eV . (39) 



7713— ~ V™si + l^l 2 (Amjj + Att^ 2 ) ~ 2.57 (1.39) eV . (40) 



(777) 



cos 2 6» jAm 2 sX + sin 2 y Atti^ e ic * 2 



Atti 2 ! a/1 - sin 2 20 sin 2 a 2 /2 . 



(41) 



All these expression are almost identical to the ones for scenarios NSS and ISS, because 
the leading contributions to all observables correspond to a situation with three quasi- 



degenerate active neutrinos having a mass yAm^. Therefore mass-related observables 
can not distinguish these cases, unless the mass-squared differences Am s 2 j and Atti^ are 
very much different from each other. 

6.4 Scheme SISb 

In this scheme we have 7774 > 7772 > mi > 7773 > 7775, i.e., the sterile neutrinos are above 
and below three inversely ordered active neutrinos, see Fig. [31 One finds 

Att^ = m\ — m\ and Am 2 2 = m\ — m\ , (42) 



14 



and can express the individual masses as 



m3 = J Am 2 , + mj , rrii — \ ^ m 1i + Am? + m| , 

, z (43) 

m 2 = J Am s \ + Am 2 A + Am^ + m§ , m 4 = w Am^ + Am^ + m| . 

We do not give the expressions for E, mp or (m), because the results are indistinguishable 
from scenario SNSb and therefore also from NSS and ISS. Again, matter effects in neutrino 
oscillation experiments could be used to distinguish between the scenarios. 



7 Discussions and Summary 

Adding two sterile neutrinos to the three active ones gives rise to eight possible mass 
orderings, out of which the right one should be identified in order to pin down the flavor 
structure of the neutrino mass matrix. We have investigated how and if mass-related 
measurements can do the job. In addition, we studied the general properties of the non- 
oscillation observables in scenarios with two sterile neutrinos. The possible mass orderings 
are shown schematically in Figs. [fl [2] and [31 Apart from the usual 3-generation masses 
and mixing parameters we have to cope with two additional mixing matrix elements \U e ^\ 
and \U e ^\ as well as with two mass-squared differences Am s \ and Am s 2 2 . Without loss of 
generality we can assume Am^ > Am 2 2 and associate Am s \ with the state 5 and Am 2 2 
with state 4, respectively. 

We use the following nomenclature for the eight different schemes: 

(i) SSX, where X = N for a normal and X = I for an inverted ordering of the mostly 
active neutrinos. In these schemes the two predominantly sterile neutrinos are heavier 
than the three predominantly active neutrinos (2+3 scenarios); 

(ii) XSS (X = N or I as before), where the two predominantly sterile neutrinos are lighter 
than the three predominantly active neutrinos (3+2 scenarios); 

(iii) SXS with X = N or I, where the three active neutrinos are sandwiched between the 
sterile ones (1+3+1 scenarios). In this class there can be four possible scenarios 
which we denote as SXSa and SXSb. The scheme SXSa corresponds to the state 5 
higher than the three active states and SXSb corresponds to the state 5 lower than 
the three active states. Those scenarios are strictly speaking not covered by the 
available analyzes of scenarios with two sterile neutrinos. In absence of any fit of this 
possibility, we assumed for simplicity that the parameters are the same as for the 
other scenarios. 

The following general comments can be made about the different mass related observable^: 
The sum of neutrino masses depends basically only on the new mass-squared differences, 

Predictions for mass-related observables in the presence of yet another sterile neutrino are discussed 
in Appendix [XI 
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mp (eV) 


(m) (eV) 


E(eV) 


SSN 


0.15 -T- 0.52 


0.0 4- 0.11 


3.33^ 3.85 


SSI 


0.16 -T- 0.52 


0.0 4 0.16 


3.33 4 3.85 


NSS, ISS, SNSb, SISb 


2.3 2.7 


0.43 2.72 


9.16 -T- 10.80 


SNSa, SISa 


0.89 4 1.11 


0.09 4 1.03 


5.19 4 5.92 


Current Bound 


2.2 eV (95% C.L.) 


~1 eV (90% C.L.) 


~1 eV 



Table 2: The ranges of the predictions for mp, (m) and £ according to the various mass 
orderings for m smallest < 0.1 eV and assuming that Am s 2 j = Am^ and Am^ = Arh^ 2 . Also 
shown are the current bounds on these observables. 



and typical (minimal) values are J Am^ + yAm^, 3^/Am^ or 3yAmg 2 , depending on 
the mass ordering. Given the best-fit values and allowed ranges of masses this is already 
in conflict with the standard cosmological scenario, as discussed below. 
The parameters relevant for neutrinoless double beta decay and direct beta-decay searches 
can be written as a contribution from the three mostly active states and the two mostly 
sterile neutrinos: 

(m) = (m) 3 + (m) st and mp = J (rnffi + (m s ^) 2 . (44) 

where |(m) 3 | and are the expressions known from 3-flavor analyzes, see Table [TJ All 
cases are possible in Eq. ( 1441) : dominance of the sterile contribution, dominance of the 
active contribution, and equal-sized contributions, leading (only in (m)) potentially to 
complete cancellation. 

In general, the mass-related observables can not distinguish between a normal or inverted 
ordering of the three active neutrinos, with the exception of schemes SSN and SSI, which 
have different predictions for (m). It will however be difficult to test this difference in 
practise, as precise knowledge of the oscillations parameters is required. Only in the 
schemes SSN and SSI the magnitude of \U ei \ and \U e5 \ is crucial for the predictions of 
(m) and mp. For all other mass orderings the dependence on |{7 e4 | and \U e §\ is suppressed 
(S does not depend on \ U e ^\ and |£/ e 5|)- Scenarios SNSb and SISb are indistinguishable from 
scenarios NSS and ISS if the mass-squared differences are equal or very similar. It turns out 
that in order to summarize all phenomenology of the mass-related observables it suffices 
to plot them for four schemes: SSN, SSI, SNSa (covering also SISa) and NSS (covering 
also ISS, SNSb, SISb). Interestingly, these four cases have also the same phenomenology 
in what regards decays of high energy astrophysical neutrinos, see Appendix [B] In Table [2] 
we present for the parameter ranges given in Eqs. ([3]) and (jlj) the predictions of the three 
quantities mp, (m), S for the four types of mass orderings in the realistic case when the 
smallest neutrino mass is smaller than 0.1 eV. Also shown are the current bounds on these 
observables. Table [3] summarizes what the various schemes mean for KATRIN and for 
future experiments searching for neutrinoless double beta decay. The following conclusions 
can be drawn: 
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scheme 


feature 


KATRIN 


Ov/3/3 


SSN 


NH plus u si , v S2 


maybe 


maybe 


SSI 


IH plus u Sl , v S2 


maybe 


maybe 


NSS, ISS 


QD with \/Am 2 sl 


yes 


yes 


SNSb, SISb 


QD with \jAm 2 sl 


yes 


yes 


SNSa, SISa 


QD with sjAm 2 2 


yes 


yes 



Table 3: The various schemes with two sterile neutrinos and their meaning for KATRIN 
and future Qv(3/3 experiments. We assumed that Am^ and Am 2 2 are larger than 0.1 eV 2 . 



• scenarios SSN and SSI predict for all observables the smallest values. Hence, they 
are the easiest to rule out. The other scenarios correspond at leading order to quasi- 
degenerate 3-neutrino scenarios with the common mass scale given by ^J~Krn^ 1 (NSS, 

ISS), \/Am 2 s2 (SNSa, SISa) or ^Am 2 sl (SNSb, SISb); 

• model independent constraints on neutrino masses stem from direct searches in the 
spectra of beta-decays. The Mainz data give the constraint nip < 2.2 eV at 95% 
C.L. (Ax 2 = 4) [2T]. Following the procedure in [3T] the 68% C.L. (A^ 2 = 1) bound 
is 0.7 eV and the 99.73% C.L. (Ax 2 = 9) bound is 4.0 eV. Scenarios SSN and SSI 
with both sterile neutrinos being heavier than the active ones can have unobservably 
small mp and are consistent with the current bound even at la as can be seen from 
Table 2. The Table also shows that scenarios SNSa and SISa are allowed at 2cr, 
while the scenarios NSS/ISS/SNSb/SISb are consistent with the Mainz result at 3a. 
The other six mass orderings will definitely result in a signal in KATRIN. If the two 
sterile neutrinos are lighter than the active ones (NSS and ISS), then there is a direct 
correspondence nip ~ <J Am^ , which can be used to rule out part of the allowed range 

of Am 2 x already at the current stage. Scenarios SNSa and SISa predict nip ~ Afh 2 2 , 
which will rule out part of its allowed range in the near future; 

• the effective mass governing neutrinoless double beta decay can only vanish when the 
sterile neutrinos are heavier than the active ones (SSN and SSI, for which the largest 
value is 0.1 eV). In the other orderings the non-maximality of the solar neutrino 
mixing angle renders (m) non-zero. If the sterile neutrinos are lighter than the active 
ones (NSS and ISS), then (m) is larger than 0.2 eV and has a maximal value above 
the current limit^of ~ 1 eV at 90% C.L. Consequently, this scenario can be ruled out 
by a stronger limit on (m) and one can also constrain parameters with the current 
limit. This concerns in particular sin 2 6@ and the Majorana phase a<i- Unfortunately, 
the sterile neutrino parameters are such that this scenario is hardly distinguishable 

8 Note that in our analysis we are not using the data coming from the positive evidence claimed by a 
part of the Heidelberg-Moscow collaboration [32] and thus we have only an upper bound. 
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from the scenarios in which one sterile neutrino is heavier and the other one lighter 
than the active ones. Whereas telling apart NSS/ISS from SNSb/SISb is basically 
impossible if Am^ is similar to Am s 2 b distinguishing SNSa/SISa from NSS/ISS re- 
quires a condition of the form JAm^ cos 29 Q > £ J Am^, where ( denotes the 
nuclear matrix element uncertainty. 

In general, the sum of neutrino masses in the scenarios under study is always larger than 
about 2 eV. This minimal value is obtained in the two schemes SSN and SSI in which 
the sterile neutrinos are heavier than the active ones. In addition, the two mass-squared 
differences related to the MiniBooNE/LSND experiments should be rather small, because 
S = \J Am s 2 ! + \J Am^- The other scenarios have sizable values of S, approaching up to 10 
eV if the sterile neutrinos are lighter than the active ones. However, cosmological limits 
can always be evaded or relaxed. Nevertheless, we add some discussion on typical limits 
obtained in the literature (for an overview, see p2]), which typically constrain both the 
sum of neutrino masses and the effective number of neutrino species iV e fj- contributing to 
the radiation density. Those limits depend also on N m , which is the number of equally 
massive species. As one example, we focus on Ref. [T3], in which likelihood contours are 
provided in the iV e fj-E plane for three cases: (i) N m = iV eff , (ii) N m = 3 and (iii) N m = 1. 
The value N e g = 5 is allowed only at about 99% C.L. in all the above cases and the bounds 
on S are 0.62, 0.57 and 0.41 eV, respectively (all at 95% C.L.). The type la supernova data 
from SNLS, large scale structure data from 2DF and SDSS, baryon acoustic oscillation data 
from SDSS, CMB anisotropy data from WMAP and the smaller scale measurement by the 
BOOMERANG experiment were included in that analysis. Adding the Lyman-a forest 
data gives even stronger bounds [19], leaving out the baryon acoustic oscillation relaxes 
the limits [TT] . For the unrealistic case when the lightest neutrino is heavier than 1 eV 
we have N m = N e s = 5 and case (i) applies. The resulting E in our scenarios is of course 
much larger than allowed. Another example is when in scenarios NSS/ISS the smallest 
mass can be neglected. Then we have three quasi-degenerate neutrinos, but also one other 
massive neutrino with mass ^Am^ — Am s 2 2 . This possibility, as well as the other cases 
we study, is not covered by the analysis in [18] or in any other paper we are aware of. 
Nevertheless, we can safely assumed that limits on the sum of masses for N e g = 5 do not 
exceed ~ 1 eV. Consequently, and not surprisingly, all scenarios with two sterile neutrinos 
have serious problems with cosmology and require non-standard physics (primordial lepton 
asymmetries, low reheating temperature, additional neutrino interactions,. . .) as described, 
e.g., in [291 [30]. 

To conclude, scenarios with two sterile neutrinos offer rich and interesting phenomenology. 
In six of the eight allowed cases KATRIN and future Qv/3/3 experiments will find a signal. 
Mass-related observables alone, however, can not identify the correct mass ordering com- 

9 A more quantitative and accurate estimate on the joint constraint on the number of neutrino species 
and the sum of neutrino masses from cosmology in the various scenarios and estimating the A% 2 would 
require a more detailed analysis of the cosmological data sets which is clearly not in the purview of the 
present analysis. 
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pletely, which leaves room for further studies in order to disentangle the possibilities by 
means of oscillation experiments. 
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A On Scenarios with three sterile Neutrinos 

For the sake of completeness we summarize briefly the formulae and results for the mass- 
related observables in case when three sterile neutrinos are added. No qualitatively new 
aspects are found in these scenarios. The authors of Ref. j9] also performed an analysis of 
this possibility and it was found that no significant improvement of the fit can be achieved 
in this way. The best-fit values for the mass-squared differences are 

Am s 2 i = 1.84 eV , Am s 2 2 = 0.83 eV , Am 2 s3 = 0.46 eV . (Al) 

Note that there can be scenarios in which the fit of Ref. [9] does not apply. There is no 
information given on the mixing matrix elements, let us therefore take for simplicity the 
values 

\U ei \ = 0.1 for % = 4,5,6 . (A2) 

We will estimate now the values of the mass-related observables for all possible mass 
orderings. Again, we fix m li2i 3 to be responsible for the oscillations of solar and atmospheric 
neutrino oscillations. With the details given in the main text, it is quite easy to obtain 
the following formulae, which are valid when the smallest neutrino mass is neglected. As 
for two sterile neutrinos, the non-oscillation probes can not distinguish whether the active 
neutrinos are normally or inversely ordered. There are in total 16 possible mass orderings, 
which we group in 4 classes: 

(i) the first type of mass spectrum holds when the three sterile neutrinos are all heavier 
than the active ones (SSSN and SSSI). In this case, 



Am 2 sl + ^Am 2 2 + ^ Am 2 s3 ~ 2.9 eV , 
mp ~ \/\U e4 \ 2 Am 2 s3 + \U e5 \ 2 Am 2 2 + \U e6 \ 2 Am 2 sl ~ 0.18 eV 

(m) ~ 



U eA \ 2 JAm 2 3 + \U e5 \ 2 JAm 2 s2 e^ 4 + \U e& \ 2 J Am 2 sl e 



< 0.29 eV , 

(A3) 
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where 0:54 = 0:5 — and a§± = a§ — at/t are combinations of Majorana phases. The 
situation is somewhat similar to scenarios SSN and SSI; 

(ii) a second class of spectra is found when the three sterile neutrinos are lighter than 
the three active ones (scenarios NSSS and ISSS). One has 

£ ~ 3^Am s 2 ! ~ 4.1 eV , 

mp ^ \jAm 2 sl ~ 1.36 eV , (A4) 
(m) ~ \jAm 2 sl \fl - sin 2 26 Q sin 2 a 2 /2 ~ (0.54 4- 1.36) eV . 

This is similar to the possibilities NSS and ISS; 

(iii) two neutrinos can be heavier than the active ones which in turn are heavier than 
the third sterile state. There are three possibilities for this. If the two heavy sterile 
neutrinos correspond to Am^ and Am 2 2 , then the resulting schemes SSNSa and 
SSISa give 

E ~ 3 V //W[ + ^3 + Am 2 , + ^Am 2 3 + Am 2 , ~ 4.7 eV , 

~ ^3 + \U eb \ 2 (Am 2 3 + Am s 2 2 ) + \U e& \ 2 (Am 2 3 + Am 2 sl ) ~ 0.70 eV , (A5) 

(m) ~ v / Am| \A - sin2 2 #© sin2 "2/2 ~ (0.27 -=- 0.68) eV . 

Also possible is that the two heavy sterile neutrinos correspond to Am^ and Am 2 3 (schemes 
SSNSb and SSISb), in which case 

E ~ 3^/Am s 2 2 + xjAm 2 2 + Am 2 3 + ^ + A< ~ 5.5 eV , 
m^ ~ ^/Am s 2 2 + \U e5 \ 2 (Am s 2 2 + Am s 2 3 ) + \U e6 \ 2 (Am 2 s2 + Am 2 sl ) ~ 0.93 eV , (A6) 
(m) ~ V /Am s 2 2 ^1 - sin 2 29 Q sin 2 a 2 /2 ~ (0.36 -=- 0.91) eV . 

Finally, the two heavy sterile neutrinos can correspond to Am 2 2 and Am 2 3 (schemes 
SSNSc and SSISc): 

E ~ 3/Am s 2 ! + yjAm? al + Am 2 , + yj Am 2 sl + Am 2 2 ~ 7.2 eV , 
m^ ~ /Am 2 ! + |f/ e5 | 2 (Am 2 ! + Am s 2 3 ) + \U e6 \ 2 {Am 2 sl + Am 2 2 ) ~ 1.37 eV , (A7) 
(m) ~ ^Am 2 ! ^1 - sin 2 29 Q sin 2 a 2 /2 ~ (0.54 4 1.36) eV . 

(iv) the fourth class of mass orderings are scenarios in which the two sterile neutrinos are 
lighter than the active ones which in turn are lighter than the last sterile one. As 
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for the previous class of scenarios, three possibilities are present. In scenarios SNSSa 
and SISSa the two light sterile neutrinos correspond to Am 2 2 and Am s 2 3 : 

E ~ 3^Am s 2 2 + ^Am s 2 2 - Am s 2 3 + y/Am 2 al + Am 2 s2 ~ 5.0 eV , 
mp ~ ^Am a 2 2 + |t/ e5 | 2 (Am s 2 2 - Am s 2 3 ) + |C/ e6 | 2 (Am 2 sl + Am 2 2 ) ~ 0.93 eV , (A8) 
(m) ~ y / Am| /l - sin 2 2# sin 2 a 2 /2 ~ (0.36 -r- 0.91) eV . 

If the neutrinos associated with Am s \ and Am 2 3 are lighter (scenarios SNSSb and 
SISSb), then 

S ~ 3^/Am 2 ! + ^Am s \ - Am 2 3 + ^Am s \ + Am 2 2 ~ 6.9 eV , 
m/3 ~ /Am 2 ! + |t/ e5 | 2 (Am 2 ! - Am 2 3 ) + \U e6 \ 2 (Am 2 sl + Am 2 s2 ) ~ 1.37 eV , (A9) 
(m) ~ A/Am 2 ! \/l — sin 2# sin 2 a 2 /2 ~ (0.54 -J- 1.36) eV . 

Finally, in scenarios SNSSc and SISSc the heaviest neutrino corresponds to Am 2 3 : 

S ~ 3^/Am 2 ! + ^Am s 2 i - Am 2 2 + ^Am^ + Am 2 3 ~ 6.6 eV , 
m^ ~ ^Am s 2 ! + |£7 e5 | 2 (Am 2 x - Am 2 2 ) + \U e6 \ 2 (Am 2 sl + Am 2 3 ) ~ 1.37 eV , 
(m) ~ /Am 2 ! yl - sin 2 2# sin 2 a 2 /2 ~ (0.54 1.36) eV . 

(A10) 

Hence, except for scenarios SSSN and SSSI one can expect a signal in KATRIN and in 
neutrinoless double beta decay searches. The latter has in this case only three main pre- 
dictions, given by a quasi-degenerate scenario with a common mass scale ^ Am 2 x , Am 2 2 

or J Am 2 3 . In order to distinguish these cases, conditions in analogy to Eq. flM|) should be 
fulfilled. Note that strictly speaking the analysis of Ref. [9j does not apply to the scenarios 
discussed in items (iii), (iv) and (v). The reason is the same as the one discussed in Section 
m for schemes SNSa/SISa/SNSb/SISb. 

B An alternative non-oscillation probe of Neutrino 
Spectra: Decay of astrophysical Neutrinos 

Another interesting non-oscillation probe to distinguish different mass orderings is the 
decay of astrophysical neutrinos [33]. Albeit such an analysis depends crucially on the 
non-trivial assumption that neutrinos decay, it is an interesting exercise to investigate the 
implications of such a situational. Different astrophysical sources can generate a neutrino 

10 For other uses of sterile neutrinos in neutrino telescopes, see [34] . 
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flux with a certain initial composition in u e , and v T . Assuming that all neutrino states 
except for the lightest Ui decay, leads to 



$ e : : $ T = \U ei \ 2 : \U^\ 2 : |£/ Tl | 2 , (Bl) 

where <3> Q with a = e, /j,, r is the flux of neutrinos and anti-neutrinos of flavor a which 
reaches Earth. Note that in the decay scenario this final flux is independent on the initial 
flavor composition. The crucial observation is that the surviving lightest neutrino state 
being z/j, the flux $ Q , which is proportional to |t/m| 2 , may differ in the eight neutrino mass 
orderings under study. For scenario SSN we have i — 1, while for SSI it holds % = 3. This 
corresponds to the situation in the three-flavor scenarios studied in Refs. [33J. In contrast, 
for scenarios NSS, ISS, SNSb and SISb we have i — 5, while in the orderings SNSa and 
SISa it holds that % = 4. There are therefore four different possibilities, and the mass 
orderings sharing the same phenomenology are in fact the same as the ones sharing the 
same phenomenology of the mass-related observables. 

What is eventually measured in neutrino telescopes like IceCube [35] are ratios of fluxes, 
and here for illustrative purposes we will focus on the ratio 



Rep = I 1 , (B2) 



which can be obtained by comparing the rate of shower and muon events [36J. Taking 
ratios including v T into account will complicate the situation considerably, as the mixing 
elements |Z7 T 4| and \U T ^\ enter the game, which are basically unconstrained (one could 
however use these ratios to obtain information on these elements). We assume maximal 
atmospheric neutrino mixing, take the best-fit values from Eqs. ([31 HD and use from Ref. [9] 
that \U^\ = 0.12 as well as \U^\ = 0.16. It follows for the ratios that i?^ N = 2/ tan 2 9 Q = 

4.7, R*f = 2sin 2 £ C HOOZ = 0, R™ S > ISS, SNSb, SISb = \ Ue tf/\U^ = 1 and ^SNS^SISa = 

I 1 V I 1 2 — 0.47. The ratios are easily distinguishable from each other, in particular 
SSN and SSI, which we have shown in Section H] to be very similar in the mass-related 
observables. Note however that for standard astrophysical sources (an initial composition 
of 1 : 2 : 0) and no decay the ratio R efl is equal to 1 for maximal atmospheric mixing 
and U e 3 = 0, i.e., identical to _R^f s ' ' b ' b . In addition, taking the uncertainty of 
the mixing matrix elements into account complicates the situation further. Taking the 3a 
ranges from Ref. [16J (in particular sin 2 9\ = (0. 32-^0.64) for atmospheric neutrino mixing) 
and assuming again a 50% uncertainty on the sterile neutrino parameters gives that -R^ N = 
(2.2-10.8), R™ 1 = (0-0.13), ^nss, iss, SNSb, sisb = ( .ll-=-9.0) and #SNSa,siSa = (q.06-4.0). 
The standard scenario predicts R e ^ between 0.73 and 1.19, a range which is covered by all 
decay scenarios except for SSN and SSI. The ratios are now overlapping and except for the 
cases of measuring very small (R efl < 0.06) or large ratios (R efM > 9) no possibility can be 
unambiguously identified. However, certain cases can be ruled out, for instance scenarios 
SSI, SNSa and SISa for an observation of R efJi > 4. To fully disentangle the different cases 
the mixing parameters the mixing matrix elements should be known much more precisely. 
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Figure 1: Allowed 2+3 mass orderings which are defined by having two sterile neutrinos 
heavier than the three active neutrinos. (Not to scale). 
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Figure 2: Allowed 3+2 mass orderings which are defined by having two sterile neutrinos 
lighter than the three active neutrinos. (Not to scale). 
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Figure 3: Allowed 1+3+1 mass orderings which are defined by having one sterile neutrino 
heavier than the three active ones which in turn are heavier than the second sterile neutrino. 
(Not to scale). Note that not necessarily Am^ = Am^ and Am^ 
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Figure 4: Scenarios with three active and two sterile neutrinos: the individual neutrino 
masses as a function of the smallest neutrino mass for scenarios SSN, SSI, NSS and SNSa. 
For the mass-squared differences related to the sterile neutrinos the best-fit point given in 
Eq. (J3J) is used and we assumed that Am^ = Am^ and Am^ = Am^. Scenario ISS is 
indistinguishable from case NSS and SISa from SNSa. The schemes SNSb and SISb are 
very similar to NSS. 
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Figure 5: The sum of neutrino masses E and the kinematic neutrino mass mp for scenarios 
SSN (top left), SSI (top right), NSS (bottom left) and SNSa (bottom right). The solid lines 
give the values of the respective observable at the best-fit point Eq. (Tjjwhile the shaded 
regions are obtained by varying the parameters involved in their corresponding ranges from 
Eqs. (12]) and We assumed that Am^ = Am^ and Am^ = Am s 2 2 . Scenario ISS is 
indistinguishable from case NSS, SISa is indistinguishable from SNSa, and SNSb/SISb are 
indistinguishable from NSS. For these two observables SSN and SSI give identical results. 
Also indicated is the KATRIN sensitivity on mp of 0.3 eV. 
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Figure 6: The effective mass (m) for scenarios SSN (top left), SSI (top right), NSS (bottom 
left) and SNSa (bottom right). The magenta shaded (darker) regions correspond to the the 
best-fit point from Eq. (j4j) and allowing the phases to take arbitrary values in the interval [0 
: 27r]. The blue shaded (lighter) regions are obtained by varying the mass and mixing angles 
as well in the allowed ranges from Eqs. (ED and (jBJ). We assumed that Am^ = Am^ and 
Am s 2 2 = A?rig 2 . Scenario ISS is indistinguishable from case NSS, SISa is indistinguishable 
from SNSa, and SNSb/SISb are indistinguishable from NSS. Also indicated is the current 
limit of 1 eV and a future bound of 0.04 eV. 
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Figure 7: Same as Fig. [6] for the second typical sterile parameter point and the correspond- 
ing range from Eq. (J5J). 
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Figure 8: Allowed values of the Majorana phase a 2 and sin 2 9 Q in scenarios NSS and ISS 
for the current limit on (m) of 1 eV (left) and a future limit of 0.5 eV (right). The 3a 
ranges from Eq. ([3]) and the ranges around the best-fit point from Eq. (J3J) are used. 
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